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RSTN  Normalized  displacement  response  from  Bredtng  (1982)  s 

LP:  long  period  band  recording 

MP:  middle  period  band  recording 

SPK:  short  period  band  KS-16000  recording 

SP7  :  short  period  band  S-750  recording 

Location  of  events  examined  in  this  report  (see  Table  T)  1 

and  RSTN  station  locations  Mn  and  M  are  shown  for  each 

n 

event . 

NTS-to-RSSD  (11.49°)  from  top  to  bottom:  Short  period  ^ 

Z,N,E  (40  seconds)  Mid  period  Z,N,F  (175  seconds)  and  Long 
period  Z,N,E  (1500  seconds).  The  vertical  scale  tick 
marks  are  lOOnm  for  each  trace,  except  the  long  period  Z. 

The  tick  marks  for  the  LPZ  are  1  micron  increments.  The  P 
wave  window  for  spectral  estimation  In  Figure  4  Is  shown 
on  each  seismogram  Indicated  by  the  vertical  bars  labeled 
"start”  and  "end".  Noise  windows  were  selected  prior  to 
the  P  arrival.  The  P  character  of  P  Is  clear  on  the 
short  period  Z  with  strong  radial  polarization.  Lg  Is 
clearly  observable  on  the  Mid  period  hand  and  I.R  Is  well 
recorded  on  the  long  period  band.  The  P  wave  Is  not 
detectable  on  the  long  period. 

NTS-to-RSSD  (11.49°)  short  period  1-component  P  wave  8 

amplitude  spectra  with  noise  spectra  superimposed.  Units 
are  nm-(sec)  at  1  Hz.  Spectra  are  not  corrected  for 
frequency  dependent  Instrument  response.  Slgnal-to-nolse 
ratio  1 s  estimated  as  450,  250,  and  525  on  Z , N , F 
respectively  all  at  1.1  Hz. 

NTS-to-RSSD  (11.59°)  Mid  period  1-component  amplitude  9 

spectra  with  noise  spectra  superimposed.  Slgnal-to-nolse 
1 8  estimated  as  566,  150,  and  150  on  Z,N,E  at  .98,  1.12, 
and  .9  Hz  respectively.  Spectra  are  asymptotic  to  noise 
levels  at  frequencies  below  0.20  Hz  ind^^tlve  of  a  long 
period  P  wave  null.  Units  are  nm-(sec)  at  1  Hz. 

Spectra  are  uncorrected  for  frequency  dependent  instrument 
response . 

Calone-to-RSSD  (14.78°)  (Scales  same  as  Figure  1)  P  wave,  11 
Love  and  Rayleigh  surface  waves  are  the  long  period 
records . 

Calone-to-RSSD  (14.78°)  Mid  period  P  wave  amplitude  1? 

spectra  with  noise  estimates.  P  wave  spectra  remain  above 
the  noise  across  the  mtcroselsm  band  at  0.2  Hz. 

Calone-to-RSSD  (14.78°)  Long  period  P  wave  amplitude  ] 1 

spectra  with  noise  estimates.  P  wave  spectra  remain  above 
the  noise  at  frequencies  below  0.1  Hz. 
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NTS-to-RSSD  (11.49*)  Long  period  L_  amplitude  spectra  with 
noise. 
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period  and  mid  period  Z&E.  (Azimuth  of  approach  ESE)  S 
waves  and  Lg  are  well  developed  on  the  mid  period  band.  A 
Rayleigh  pulse  is  prominent  on  the  long  period  band. 

Dakota-to-RSSD  (3.3*,  4.4)  Short  period  P  wave 

amplitude  spectra  with  noise  levels.  P  wave  spectra  above 
the  noise  throughout  the  band  width.  Signal  to  noise  has 
multiple  local  maxima  at  3.7  and  13.1  Hz.  (Noise  estimate 
la  plotted  1  decade  too  high) 
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spectra  with  maximum  slgnal-to-nolse  near  3  Hz. 

NTS  short  and  long  period  P  wave  detection:  P  logarithmic 
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values  Is  for  short  period  P  and  the  bottom  pair  of  values 
is  for  long  period  P.  The  M,  5.7  NTS  explosion  had  a  P 
wave  logarithmic  slgnal-to-nolse  ratio  of  2.7  at  1.1  Hz 
RSSD  (3.3*).  Similarly  the  long  period  P  logarithmic 
slgnal-to-nolse  ratio  was  0.7  at  0.34  Hz.  A  detection 
threshold  of  M,  3.0  for  this  path  is  Indicated  for  P  and  M^ 
5.0  for  long  period  P. 

Dakota  earthquake  (M.  4.5)  short  period  and  long  period  P 
wave  detection.  RSNT  records  were  not  available  for  this 
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EVALUATION  OF  THE  RSTN  NETWORK 


Introduction 

The  Regional  Seismic  Test  Network  (RSTN)  is  a  5-statlon  seismic 
network  designed  to  monitor  regional  seismicity  on  the  scale  of  a  large 
continental  shield.  The  RSTN  serves  as  a  prototype  for  a  network  that 
could  be  deployed  for  test  ban  treaty  verification.  Both  regional  and 
teleselsmic  signals  are  digitally  recorded  for  three  components  of  earth 
motion  in  three  overlapping  frequency  bands  covering  0.02  to  10  Hz. 

The  instruments,  signal  conditioning,  and  recording  facilities  are 
discussed  in  Rodgers  and  Hammel  (1981),  and  Breding  (1982). 

Each  station  consists  of  a  three  component  broadband  Teledyne 
Ceotech  KS-36000  and  a  high  frequency  Teledyne  Geotech  S-750  in  a 
borehole  at  a  depth  of  100  meters.  The  KS-36000  output  signals  are 
processed  to  yield  a  long  period  band,  an  intermediate  period  band,  and 
a  short  period  band  (see  Figure  1).  The  digital  data  is  transmitted  by 
satellite  to  the  System  Control  and  Receiving  Station  in  Albequerque,  N. 
M.  (SCARS),  and  the  Center  for  Seismic  Studies  (CSS)  in  Arlington,  Va., 
as  well  as  other  sites  in  North  America.  The  data  used  in  this  report 
was  partly  received  at  CSS  and  partly  recorded  at  SCARS.  A  description 
of  RSTN  station  localities  is  contained  In  Taylor  and  Qualheim  (1983). 
All  RSTN  stations  are  located  on  Precambrian  or  Paleozoic  sediment  at 
"hard  rock"  sites  of  the  Canadian  Shield.  The  network  spans  nearly  27° 
of  latitude  and  40°  of  longitude  (Figure  2). 

Because  the  RSTN  is  a  sparse  network  of  three-component  "broadband" 
stations  it  poses  unique  problems  in  the  task  of  developing  a  regional 
system  of  automatic  detection  and  association.  The  P^,  Pg,  Sn,  Lg 
regional  phases  have  not  been  incorporated  into  a  detection,  location, 
and  association  process.  Three-component  stations  with  three  frequency 
bandwidths  offer  the  opportunity  to  determine  the  optimal  frequency 
bands  for  detection  and  identification  of  these  phases.  Regional 
detection  becomes  substantially  Important  in  the  context  of  a  CTBT 
because  of  agreements  between  the  USA  and  USSR  that  Internal  3-component 
stations  could  be  Implemented  in  the  event  of  such  a  treaty. 


Normalized  (counts/^) 


Figure  2.7  HS'l'N  Normalized  Displacement  Amplitude  Responses 


Figure  1  RSTN  Normalized  displacement  response  from  Bredlng  (1982) 
LP:  long  period  band  recording 
MP:  middle  period  band  recording 
SPKs  short  period  band  KS-36000  recording 
SP7 :  short  period  band  S-750  recording 
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Location  of  events  examined  in  thla  report  (see  Table  I)  and 
RSTN  station  locations  and  are  shown  for  each  event. 


TABLE  T 


Event  Location 


ml 

NAME 

9 

LAT(N) 

LON(W) 

OT 

DATE 

mb 

5.4 

NTS 

1 

37 

116 

14:00 

82/217 

5.7 

5.4 

CA1.0NE 

2 

16.3 

120.5 

22:26 

82/298 

4.7 

1.5 

ARK 

1 

35.1 

92.2 

02:10 

83/019 

4.1 

4.5 

YELLOW 

4 

44.5 

110.6 

20.25 

83/037 

4.2 

6.5 

CAL 

5 

36.2 

120.3 

23:42 

83/122 

5.6 

4.1 

GASPE 

6 

49.1 

67.0 

19:15 

83/017 

4.1 

2.9 

NY 

7 

41.5 

73.6 

19:59 

83/057 

1.5 

(M_ )  6.1 

MEX 

8 

16.1 

95.2 

08:17 

83/024 

5.7 

B  4.4 

DAKOTA 

9 

44.2 

99.4 

06:12 

81/065 

4.5 

TABLE  IT 

RSTN 

Station  Locations 

RSTN 

Lat 

Long 

F.lev 

RSNY 

RSCP 

RSSD 

RSON 

RSNT 

(New  York) 

(Cumberland  Plateau) 
(South  Dakota) 

(Ontario) 

(North  West  Territories) 

44 “32 ' 54 "N 
35°36'00”N 

44  *07 ' 1 3 . 5”N 
50®5l * 32"N 

62  ® 28 ' 47 "N 

74  °  3 1 '  48"W 
85°34 '08“W 

104  *02 ' 10 . 3 "W 
91*42' 08“W 
lu'is'io-w 

457m 

581m 

2060m 

405m 

191m 

The  approach  of  thla  study  was  to  select  9  events  recorded  by  the 

network  and  manually  analyze  the  recordings  using  available  software  at 

the  CSS  (Table  I,  and  Figure  2).  Various  phases  were  identified  and 

"picked’’  as  well  as  possible,  back  azimuths  were  estimated,  and 

signal-to-noiae  ratios  were  determined.  Finally,  detection  thresholds 

as  a  function  of  distance  were  estimated  for  P,  Lg  and  L_  phases  using 

K 

the  slgnal-to-noise  ratios  of  the  9  events.  Several  existing  CSS 
software  packages  were  tested  using  the  RSTN  evaluation  data..  In  the 
course  of  data  analysis  the  quality  of  RSTN  data  was  examined.  Due  to 
deficiencies  in  both  data  quality  and  available  software,  some  Intended 
tasks  could  not  be  performed  on  the  RSTN  data.  These  tasks  were  subject 
to  the  programming  and  computational  priorities  of  the  CSS  activities. 
Data  analysis  was  supported  with  programs  described  in  Lincoln  Lab 
manual  132  and  standard  UNIX  software.  Several  ’’bugs"  were  found  and 
corrected  in  the  existing  software.  Also,  short  digital  "drop-outs'' 
detected  in  the  data  required  additional  changes  to  the  existing 
software  to  locate  and  correct  the  data  drop-outs. 


Data  Analysis 


The  following  discussion  serves  to  Illustrate  the  analysis  performed 
at  the  CSS  with  examples  from  the  RSTN  data  used  In  this  report. 

Figure  3  shows  the  NTS  event,  day  217  of  1982,  recorded  at  RSSD. 

The  distance  Is  11.49*.  From  top-to-bottom  are  40  seconds  of 
three-component  short  period  data,  375  seconds  of  three~component  middle 
period  data,  and  1500  seconds  of  long  period  data.  The  P  wave  Is 
emergent  on  the  short  period  components  with  strong  radial  polarisation. 
The  measured  apparent  azimuth  from  the  horizontal  components  Is  within  2 
degrees  of  the  computed  back  azimuth  to  the  event.  Lg  is  clearly 
visible  on  the  middle  period  band  and  is  well  recorded  on  the  long 
period  band.  Although  the  P  arrival  is  well  recorded  on  the  middle 
period  band,  the  P  wave  train  is  not  detectable  In  the  time  domain  on 
the  long  period  band.  The  time  windows  used  for  spectral  estimation  of 
the  P  wave  are  Indicated  with  vertical  bars  labeled  “start"  and  “end". 
Amplitude  slgnal-to-nolse  ratios  (S/N)  were  estimated  in  the  time  domain 
and  in  the  frequency  domain.  The  results  are  tabulated  in  Appendix  I 
for  each  event,  station,  and  signal  channel.  The  time  domain  maximum 
amplitude  S/N  was  estimated  for  the  short  period  Z,  N,  and  E  components 
as  98,  73,  and  85  respectively.  The  predominant  periods  were  0.2,  0.18 
and  0.15  seconds  on  the  respective  Z,  N,  and  E  components. 

Figures  4a  and  4b  are  plots  of  the  short  period  and  middle  period 
channel  P-wave  amplitude  spectra  for  the  windows  indicated  In  Figure  3. 
The  noise  estimate  (dashed  line)  Is  computed  from  a  window  of  the  same 
length  for  each  component  directly  preceedlng  the  P  wave  arrival.  All 
spectra  In  this  report  are  shown  as  amplitude  spectra  uncorrected  for 
instrument.  The  maximum  S/N  was  estimated  from  each  plot  as  450,  250, 

I  and  525  at  1.1  Hz  on  the  short  period  band  and  466,  350,  and  350  at  .9 

to  1.1  Hz  on  the  middle  period  band.  The  long  period  band  P-wave  has  no 
significant  signal  above  the  noise.  The  P-wave  spectra  in  Figure  4b  are 
asymptotic  to  the  noise  level  at  0.2  Hz.  This  was  common  at  all  RSTN 
I  stations.  A  long  period  explosion  P  wave  null  Is  indicsted  by  the  time 

domain  and  spectral  domain  estimates  of  the  P-wsve  signal  strength. 

! 
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Figure  3  NTS-to-RSSD  (11.49  )  from  top  to  bottom:  Short  period  Z,N,E 
(40  seconds)  Mid  period  Z,N,E  (375  seconds)  and  Long  period 
Z,N,E  (1500  seconds)  The  vertical  scale  tick  marks  are  lOOnm 
for  each  trace,  except  the  long  period  Z.  The  tick  marks  for 
the  LPZ  are  1  micron  Increments.  The  P  wave  window  for 
spectral  estimation  In  Figure  4  Is  shown  on  each  seismogram 
Indicated  by  the  vertical  bars  labeled  "start"  and  "end". 
Noise  windows  were  selected  prior  to  the  P  arrival.  The  P 
character  of  P  Is  clear  on  the  short  period  Z  with  strong  n 
radial  polarization.  Lg  Is  clearly  observable  on  the  Mid 
period  band  and  LR  is  well  recorded  on  the  long  period  band. 
The  P  wave  Is  notaetectsble  on  the  long  period. 


NTS, RSSD 


For  comparison  with  the  explosion  traces  of  Figure  1 ,  the  California 
earthquake  of  10/25/82,  M  5.4  (  a,"  14.78°),  Is  shown  In  Figure  5a.  The 

w 

Initial  short  period  P-wave  appears  with  much  the  same  slgnal-to-nol se 
ratio  as  the  NTS  event,  while  the  long  period  P-wave  Is  clearly  visible 
on  the  bottom  three  traces.  Love  and  Rayleigh  waves  are  clearly  seen  on 
the  long  period  records.  S  waves  are  apparent  on  both  middle  and  long 
period  band  records.  In  Figures  5b  and  C,  P-wave  middle  period  and  long 
period  spectra  are  observed  above  the  noise  at  frequencies  below  0.1  Hz 
at  this  distance  and  magnitude  level  for  earthquakes. 

Lg  slgnal-to-nolse  ratios  are  estimated  In  much  the  same  way  an  P- 
wave  slgnal-to-nolse  estimates.  Figure  6a  demonstrates  the  short  period 
Lg  spectra  of  the  NTS  event  at  RSSD  (a  “  11.49°)  and  Figure  6b 
illustrates  the  same  on  the  middle  period  band  records.  Lg  has  similar 
amplitudes  on  all  three  components  of  motion.  The  maximum 
slgnal-to-nolse  ratio  for  Lg  was  found  at  0.5  to  0.7  Hz  on  both  short 
period  and  middle  period  channels.  Similarly,  the  long  period  L 
signals  are  presented  In  Figure  7. 

As  an  example  of  a  shorter  Intra-shield  propagation  path,  the  Dakota 
earthquake  recorded  at  RSSD  (  3.1°)  is  shown  In  Figure  8.  There  Is  a 

well  developed  Pg-Pn  separation  of  8  seconds  on  the  short  and  middle 
period  bands.  Propagation  Is  from  the  east  (87°)  and  the  P  wavetraln 
Is  best  observed  on  the  vertical  and  east  components.  S  waves  and  T.g 
are  well  developed  on  the  middle  period  band.  A  broadband  pulse  of 
superimposed  S  and  surface  waves  Is  prominent  on  the  long  period 
records. 

In  contrast  to  the  NTS  P-wave  spectra  at  RSSD  (Figure  4a)  the 
spectra  of  the  Dakota  earthquake  In  Figures  9a  and  b  are  above  the  noise 
to  the  Nyqulst  frequencies  at  20  Hz  and  2  Hz.  Some  aliasing  may  be 
present . 

At  this  point  we  should  point  out  some  aspects  of  the  short-period 
noise  spectra.  Tn  Figures  4a  and  6a  the  "noise"  spectra  are  all  nearly 
flat  between  6  and  10  Hz.  This  Is  common  for  RSSD.  Between  10  and  20 
Hz  the  level  declines  to  .01  mu.  The  spectra  have  not  been  corrected 
for  Instrument  response  and  the  high  frequency  level  corresponds  to  a 
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Figure  6B  NTS-to-RSSD  (11.  *9")  Mid  period  Lg  amplitude 
noise  estimates.  Maximum  slgnal-to-nolse  rat 
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Figure  8  Dakota-to-RSSD  (3.3*,  M  4.4)  Near  regional  event  shows  well 
developed  P  -Pg  separation  of  8  seconds  on  short  period  and 
■id  period  S&E.  (Azimuth  of  approach  BSE)  S  waves  and  Lg  are 
well  developed  on  the  mid  period  band.  A  Rayleigh  pulse  is 
prominent  on  the  long  period  band. 


SPP  DAKOTA. RSSD 


Dakota-to-RSSD  (3.3®,  M,  4.4)  Short  period  P  wave  amplitude 
apectra  with  noise  levels.  P  wave  spectra  above  the  noise 
throughout  the  band  width.  Signal  to  noise  has  multiple 
local  maxima  at  3.7  and  13.1  Hr.  (Noise  estimate  is  plotted 
1  decade  too  high) 
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Figure  9B  Dakota-to-RSSD  (3.3*,  4.4)  Mid  period  P  amplitude  spectra 

with  maximum  aignal-to-noise  near  3  Hz. 


signal  level  of  +.5  least  significant  bit  (LSB),  the  quantization  level. 
For  the  short  period  band,  this  corresponds  to  .01  m  |i  /  /llz  at  the 
reference  frequency  of  1.0  Hz.  At  highest  gain,  the  short-period  LSB  “ 
0.031  a  u  and  the  effect  of  quantization  noise  may  be  viewed  as  white 
noise  with  ras  amplitude  of  LSB/  / 12  (Openheim  and  Schafer,  1975)  or 
0.01  au  .  At  the  next  lowest  gain  level  the  quantization  noise  level  Is 
.07  au  •  For  a  larger  aaplltude  signal  that  was  predoalnately  recorded 
at  a  greater  gain  level,  the  quantization  noise  level  will  be  higher. 

Table  III  suaaarlzes  the  short-period  quantization  noise  levels  at 
the  four  gain  levels  that  are  used  with  the  automatic  gain-ranging.  To 
be  precise,  the  white  noise  level  contribution  to  the  final  spectra  Is  a 
sub  of  the  proportional  time  periods  at  each  automatic  gain  level. 

The  Halting  white  noise  level  Is  at  best  .01  m  n  at  l  Hz.  When 
corrected  for  lnatruaent  response  this  aeans  that  some  RSTN  short-period 
recordings  are  sensing  quantization  noise  and  not  local  ground  motion  at 
frequencies  above  10  Hz,  during  quiet  time  periods  (Figure  6a  noise 
levels).  Furthermore,  the  flat  signal  spectra  above  their  "noise" 
levels  at  high  frequencies  may  not  be  significant  since  the  signals  were 
recorded  at  a  higher  gain  levels  than  the  ambient  "noise"  signals. 
Consequently,  It  Is  not  clear  that  the  signals  observed  In  Figures  4a, 

6a  and  9a  are  at  or  below  ambient  noise  levels  at  high  frequencies. 
Previous  studies  of  ambient  noise  levels  at  quiet  sites  (Herrin,  1982; 
Fix,  1972)  Indicate  that  ambient  noise  levels  continue  to  decline  at 
high  frequencies. 

The  frequency  response  and  gain  levels  of  the  RSTN  Instruments  were 
designed  so  that  the  expected  ambient  ground  noise  averaged  over  1/2 
octave  bands,  was  an  rms  5  counts  up  to  10  Hz.  The  response  declines  at 
frequencies  above  10  Hz  to  ensure  that  aliasing  of  ground  and  system 
noise  Is  unimportant.  The  model  for  ambient  ground  noise  was  based  on 
known  quiet  sites  such  as  Lajltas  (Herrin,  1982).  Recent  reports  from 
NORSAR  (Bungum,  1983)  indicate  that  southeastern  Norway  has  very  similar 
noise  levels  at  10  Hz.  If  these  are  globally  representative  noise 
levels  then  it  Is  expected  that  such  quiet  sites  as  RSSD  will  show 
quantization  noise  from  time  to  time  for  frequencies  near  and  above  10 
Hz.  This  Is  not  unexpected  given  the  design  criteria  of  the  RSTN 
response. 
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Estimation  of  Detection  Thresholds 


An  effort  was  made  to  estimate  detection  thresholds  as  a  function  of 
distance  for  P,  Lg  and  L  within  the  three  frequency  bands.  Figure  10 
Illustrates  the  procedure  using  the  NTS  (M  5.4)  event  and  short  period, 

M 

P-wave  and  Lg  slgnal-to-nolae  estimates.  The  location  of  the  NTS  event 
la  shown  together  with  the  five  RSTN  stations.  At  each  station  the 
logarithm  (base  10)  of  the  maximum  slgnal-to-nolae  ratio  Is  shown  and 
the  frequency  at  which  It  was  measured.  In  the  case  of  RSNY  the  short 
period  Lg  was  not  detected.  In  the  case  of  RSCP,  the  Lg  record  was  not 
available.  The  difference  between  the  log(S/N)  and  the  local  magnitude 
provides  an  estimate  of  the  possible  detection  threshold.  For  Instance, 
an  estimate  of  M  3.0  la  derived  for  the  detection  threshold  at  RSSD  for 

u 

the  NTS  source  area  on  the  short  period  Instrument.  A  value  closer  to 
5.0  would  be  Indicated  for  short  period  Lg.  Similarly,  Figure  11 
presents  the  short  period  and  long  period  P-wave  atgnal-to-notse  ratios, 
(S/N),  for  the  Dakota  earthquake.  Figures  12  and  13  show  the  (S/N)  for 
the  Arkansas  and  New  York  earthquakes. 

The  was  used  In  these  analysis  where  available  since  It  more 
properly  measures  the  high  frequency  content  of  the  source  than  the  low 
frequency  measure  M,,  or  the  narrow  band  measure  m^.  Because  the 
velocity  response  of  the  middle  and  short-period  band  recording  Is  flat 
in  the  pass  bands  of  Interest,  the  amplitude  measurements  were  not 
corrected  to  period;  (A)  not  (A/T)  was  used  In  the  analysis. 

All  the  estimated  short-period  detection  thresholds  are  combined  In 
Figure  14  and  plotted  versus  distance.  Each  estimate  Is  coded  acordlng 
to  event  number.  If  no  short  period  P  detection  was  made  an  up  arrow 
(+)  Indicates  the  lower  limit  that  is  Implied  by  the  lack  of  detection. 

The  intra-shield  paths  of  events  (3)  Arkansas,  (4)  Yellowstone,  (6) 
Caspe,  (7)  New  York,  and  (9)  South  Dakota  may  have  systematically  lower 
thresholds.  The  non  intra-shleld  paths  of  (5)  NTS,  (2)  and  (5) 
California,  and  (8)  Mexico  could  be  expected  to  pose  additional 
attenuation  to  P-waves.  Some  of  the  Lg  (S/N)  Is  presented  in  Figures 
15,  16,  17,  and  18  together  with  the  Lg  (S/N)  data.  The  Lg  detection 


Arkansas  earthquake  (M  3.5)  Short  and  long  period  P  vave 
detection  and  logarlthnlc  elgnal-to-nolse  ratios.  The  event 
was  below  detection  threshold  at  RSNT  for  both  short  and  long 
period  P  waves  and  below  long  period  threshold  for  all 
stations.  Short  period  detection  threshold  of  3.8,  3.7  and 
3.3  are  Indicated  for  the  path  to  RSCP  (5. A*),  RSNY  (16*)  and 
RSON  (16*)  paths. 


Figure  19  Lg  detection  threshold  versus  distance  for  the  nine  studied 
events.  Lg  detection  thresholds  are  not  significantly  below 
P  wave  detection  thresholds  of  Figure  14. 


estimates  are  gathered  In  Figure  19  as  a  function  of  distance.  The 
estimated  Lg  detection  thresholds  are  significantly  better  than  the 
estimated  P-wave  detection  thresholds  of  Figure  14  for  distances  greater 
than  15  degrees.  This  may  be  an  artifact  that  events  from  the  tectonic 
provinces  of  North  America  dominate  the  longer  paths  sampled,  and  the 
larger  magnitude  events. 


L  detection  thresholds  derived  from  (S/N)  estimates  are  compiled  In 
R  -1  /2 

Figure  20.  The  r  curve  Is  plotted  for  comparison.  The  large  (m. 
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6.3)  Mexican  earthquake  may  yield  a  biased  estimate  of  the  L  detection 

threshold  because  of  source  scaling  effects.  Larger  events,  with  corner 

frequencies  below  1  Hz,  are  not  adequately  characterized  by  the  1  Hz 

estimate.  Since  the  m.  estimate  for  such  an  event  Is  probably  an  under 

b 

estimate  of  the  low  frequency  L  generated  by  the  event,  the  L 

K  H 

thresholds  for  (8),  the  Mexican  earthquake,  are  biased  low. 


consequently  biased  estimates  of  detection  threshold. 


Data  Flow  at  the  CSS 


Tapes  of  the  RSTN  data  In  SCARS  format  were  directly  recorded  at  the 
CSS  or  were  obtained  from  SCARS  by  request*  It  was  found  that  existing 
software  at  the  CSS  did  not  properly  Identify  and  process  calibration 
sequences.  Furthermore,  dead  spots  were  found  in  the  digital  data 
recorded  at  the  CSS.  These  data  drop-outs  were  usually  brief  intervals 
of  zeros  (c.f.  Figures  21a  and  21b).  To  alleviate  the  problem  of 
drop-outs  and  calibration  sequences,  some  tapes  were  requested  from 
SCARS  to  supplement  the  locally  recorded  data.  The  requested  SCARS 
tapes  had  significantly  fewer  digital  drop-outs,  but  the  problem  still 
persisted.  The  data  drop-outs  appear  to  be  the  result  of  variations  In 
the  satellite  signal  characteristics.  A  detection  algorithm  for  these 
brief  strings  of  zeros  was  added  to  the ‘ Transform  program  resident  at 
the  CSS. 

The  RSTN  data  was  passed  through  two  programs  NSSREAD  and  NSSTODB  to 
the  INGRES  data  base  management  system.  A  shell  script  was  developed  to 
extract  the  RSTN  digital  data  and  associated  calibration  Information  to 
a  local  data  base  consisting  of  seismograms  and  directory  information. 
Programs  display,  transform  and  spectrum  could  then  be  used  to  plot  the 
waveforms,  and  calculate  Fourier  spectra  of  the  digital  waveforms  on  a 
Tektronix  storage  tube  terminal.  The  plotting  format  options  of 
spectrum  were  modified  to  accomadate  the  desired  spectral  plots.  A 
serious  "bug"  in  the  program  transform  was  discovered  that  resulted  In 
severe  aliasing  of  the  resultant  spectra.  Corrections  were  made  to  the 
program  transform  and  tapering  options  were  added  to  the  algorithm.  The 
post  detection  arrival  detector  was  found  to  be  Inadequate  to  the  task 
of  automatic  detection.  The  algorithm  stops  with  its  first  detection 
and  could  not  be  modified.  Data  drop-outs  and  noise  bursts  confused  th<^ 
algorithm,  and  produced  too  many  false  alarms.  The  algorithm  could  not 
adjust  1 1 8  detection  criteria  to  a  changing  noise  environment. 

Amplitudes  and  predominant  periods  of  the  seismic  phases  were  measured 
with  the  disp  program  on  the  Seismic  Analysis  Station  (SAS). 
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Figure  2  2  Data  Flow  at  the  CS1! 
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FURTHER  IMPROVEMENTS  TO  AUTOMATIC  ASSOCIATION  (AA) 


Using  the  artificial  data  we  had  created  for  the  International  Data 
Exchange  Exercise,  we  tested  Automatic  Association  (AA)  with  the 
objectives  of  comparing  Its  performance  to  Center  for  Seismic  Studies 
(CSS)  AA,  and  to  develop  and  Include  some  AA  features  used  by  the 
Swedish  AA  and  the  CSS  AA  not  yet  present  In  the  Seismic  Research 
Center's  (SRC)  AA.  In  addition,  we  developed  and  Implemented  the 
following  features  Into  SRC's  AA. 

.  Kinematic/dynamic  use  of  amplitude  data 
.  two  array  event  location 

.  use  of  pP  and  sP  to  attempt  to  automatically  constrain  depth 
.  slowness  check  on  all  phases  (not  just  P  and  PKP) 

.  show  unassociated  arrivals  In  final  summary  bulletin 
.  Include  confidence  ellipses  In  summary  bulletin 

A  comparison  of  the  two  AA's  from  SRC  and  CSS  has  not  been 
accomplished  yet,  mainly  because  the  two  AA  programs  were  executing  on 
different  computers  and  a  common  mutually-readable  data  set  has  not  been 
defined. 

Documentation  for  AA  was  completed  and  printed  in  Seismic  Research 
Information  System  (SRIS),  (1982). 


Klnemat lc/Dynamic  Use  of  Amplitude 


Theoretical  basis  for  this  work  has  been  described  In  a  memorandum 
by  R.  R.  Rlandford  dated  7  October  1982  and  a  practical  embodiment  of 
the  scheme  Is  described  In  another  memorandum  by  Blandford  dated  19 
April  1983.  (Both  memoranda  are  contained  In  Appendix  TT).  An 
efficient  Newton-Raphson  procedure  for  computing  the  maximum  likelihood 
magnitude  (Ringdal,  1976 )  offers  a  computational  Improvement. 

For  a  chosen  network  of  stations,  not  necessarily  the  full  network 
of  detecting  station,  a  maximum  likelihood  estimate  of  body  wave 
magnitude  m^  using  the  amplitudes  of  the  locating  arrivals.  Clearly 
there  are  three  special  cases  to  consider  regarding  the  network  of 
"reliable”  stations,  the  stations  that  actually  detected,  and  the 
amplitudes  of  the  detections. 

1)  A  station  does  not  report  a  detection  even  though  the  m. 

n 

Indicates  that  It  should  have  detected.  2)  A  station  reports  a 
detection  and  the  m^  Indicates  that  It  should  not  have  detected.  3)  A 
station  reports  a  detection  but  the  amplitude  Is  either  too  high  or  too 
low  according  to  the  m^.  Figure  23  shows  the  the  numerical  modification 
to  the  accumlated  screen  points  for  these  three  cases,  and  Figure  24  Is 
an  up  to  date  list  of  the  numerical  values  for  the  attributes  considered 
In  the  SCREEN  subroutine. 
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•  SELECT  "RELIABLE  STATIONS",  GET  DETECTION  THRESHOLDS 


•  COMPUTE  RINGDAL  m„ 


•  IF  A  STATION  DID  NOT  DETECT  BUT  SHOULD  HAVE, 
SUBTRACT  POINTS,  i.e.: 

IF  p  (det)  *  .95:  -1,0 

p  (det)  >  .99:  -2,0 


•  IF  A  STATION  DID  DETECT  AND  REPORTS  AN  'UNUSUAL' 

AMPLITUDE,  Aq,  DISCARD  THAT  AMPLITUDE  FOR  THIS  EVENT, 


IF  p  (amp  >  Ao>  <  .01  OR 

p  (amp  <  A0)  <  .01  THEN  REJECT  A0,  START  SCREEN  OVER 


•  IF  A  STATION  DID  DETECT  BUT  SHOULD  NOT  HAVE,  DISCARD 
ALL  POINTS  DUE  TO  THAT  STATION,  i.e.: 

IF  p (det )  <  .01,  -  (ALL  POINTS  DUE  TO  STATION) 


Figure  23  Amplitude  Modifications  to  SCREEN 
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10  POINTS  IN  A-SEISMIC  ZONE 

Figure  2 A  SCREEN  Points 
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Depth  Constraint  with  pP  and  sP 


Our  experience  with  analyst-named  phases  indicates  a  tendency  to 
name  the  teleseiamic  phase  next  arriving  after  P  as  pP.  In  fact,  the 
next  arriving  phase  may  not  be  pP ,  but  it  may  be  sP,  PcP  or  some 
unidentified  phase.  Restraining  depth  to  the  analyst's  P-pP  time  during 
location  convergence  could  lead  to  difficulties  in  correct  event 
location  and  association  of  phases. 

These  considerations  lead  to  the  "soft”  treatment  of  pP  and  sP  named 
phases  described  here. 

The  list  of  P  arrivals  that  is  about  to  be  input  to  the  HYPO 
location  program  is  examined  in  time  order.  The  arrival  queue  is 
searched  for  a  detection  at  the  same  station  that  is  labeled  pP  or  sP. 

If  an  sP  1 8  found  the  search  terminates  and  the  depth  for  P-sP  is  used 
as  an  input  depth  to  HYPO.  HYPO  is  not  constrained  to  this  depth , 
however.  If  no  sP  is  found,  the  smallest  P-pP  time,  if  any  are  found, 
is  used  to  establish  the  input  depth  to  HYPO,  again  unrestrained  to 
depth.  The  Idea  here  is  to  try  to  Improve  the  initial  depth  estimate  of 
the  HYPO  solution. 

The  situation  where  the  depth  is  actually  constained  by  labeled  pP 
or  sP  phases  occurs  as  follows.  When  an  event  has  been  made  and  fully 
associated,  its  detections  are  examined  to  see  if  there  is  an  sP  phase, 
labeled  as  such  by  AA,  which  had  also  been  labeled  sP  by  the  analyst. 

If  one  is  found  this  P-sP  depth  is  used  to  constrain  HYPO  and  the  event 
is  run  again.  If  no  sP  is  found,  the  same  algorithm  is  applied  to  P-pP 
depth,  using  the  smallest  P-pP  time,  if  any.  This  algorithm  works  very 
well  with  the  artificial  data  because  analyst  labeled  "sP”  is  always 
true  in  this  data  set.  More  work  has  to  be  done  on  this  algorithm  to 
expand  its  usefulness  to  real  data. 


Other  AA  Improvements 


In  addition  to  the  two  major  Improvements  described  above  some 
continuing  enhancements  to  AA  were  accomplished*  The  slowness  check  has 
been  extended  to  include  all  phases  recognized  by  AA.  Previously  the 
slowness  check  was  limited  to  P  and  PKP  phases  during  the  location 
phase.  Since  the  travel  times  for  the  remslning  phases  are  computed  by 
polynomial  evaluation  rather  than  by  table  lookup,  this  necessitated  the 
creation  of  a  complete,  new  program  component  in  which  there  is 
duplication  of  several  of  the  computational  modules  used  in  travel  time 
determinations  during  location. 

Two  array  station  trial  epicenters  have  been  added,  based  on  the 
intersection  of  their  respective  confidence  ellipses.  SRC's  AA  uses 
this  as  the  last  pass  through  the  signal  queue,  so  in  general,  very  few 
events  are  started  this  way.  It  has  been  added  to  be  compatible  with 
other  AA  programs.  Equations  for  the  two  array  computation  are  shown  in 
Figure  25. 

Confidence  ellipse  information  has  been  added  to  the  summary 
bulletin  in  accordance  with  standard  practice. 

The  format  of  the  summary  bulletin  has  been  expanded  so  that 
"left-over"  or  unassociated  detections  are  printed  in-line  with  the 
events  and  their  association.  This  printing  is  a  great  help  in 
reviewing  the  bulletin  and  the  performance  of  AA.  Various  printing 
highlights  are  used  to  make  the  unassociated  detections  differentiate 
from  the  associated  detections. 
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Figure  25  Two-Array  Location 


Future  Improvements 


As  mentioned  earlier,  further  work  needs  to  be  done  using  pP  as  a 
depth  constraint  with  real  data.  Methods  have  to  be  worked  out  for 
making  best  use  of  the  depth  phase  Information  In  the  signal  arrival 
queue,  and  probably  for  making  a  choice  of  different  algorithms.  In 
some  cases  no  uses  at  all  would  be  made  of  pP  and  sP,  and  automatic 
methods  to  determine  that  this  would  be  the  best  strategy  needed. 

We  believe  that  the  capabilities  of  AA  ought  to  be  expanded  to 
include  and  make  use  of  regional  phases.  There  are  six  local  phases: 

Pn,  Pg,  P,  Sn,  Sg  and  S.  The  program  ACLOC  uses  regional  phases  as  well 
as  array  Information  and  other  later  phases  for  location.  ACLOC  would 
replace  HYPO  in  AA.  Conceptually,  this  is  fairly  straight-forward,  but 
considerable  programming  effort  is  required  to  make  the  resulting  large 
package  manageable. 
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MEMORANDUM 


TO: 

A.  Kerr,  W.  Dean, 

J.  Goncz,  R. 

North,  R.  Slunga,  R.  Shumway 

FROM: 

R.  Blandford 

SUBJECT: 

Dynamic-Kinematic 

Criteria  for 

Event  Reality 

DATE: 

7  October  1982 

We  now  have  a  complete  "system"  for  dynamic-kinematic  amplitude 
checking;  attached  to  this  memo.  I  think  this  system  could  be  checked 
out  fairly  well  with  the  artificial  data  John  Goncz  and  I  are 
generating.  The  subroutine  amp  may  diverge  if  very  large  or  small  TD  “ 

1  type  values  are  in  the  data.  Bob  Shumway  and  I  will  soon  distribute 
an  improved  version  of  amp.  I  think  this  code  is  in  a  "final"  state 
suitable  for  use  in  existing  AA  programs.  Only  the  I/O  needs  to  be 
modified,  and  this  is  fairly  well  isolated  in  the  event. f  and  mkpkin.f 
files . 

As  long  as  we  are  talking  about  AA,  I  propose  a  measure  of  AA 
quality,  "Association  Percentage"  defined  as  follows.  Define  events 
which  are  to  be  detected;  e.g.  events  with  5  non-array  stations;  plus 
events  with  two  arrays  and  one  confirming  station,  plus  events  with  2  P 
waves  and  one  S,  etc.  All  arrivals  in  this  set  of  events  when  perfectly 
analyzed  constitute  the  "base  set  of  arrivals".  Now  we  go  down  the 
perfect  bulletin  and  pick  out  the  first  event.  Is  there  an  event  in  the 
AA  with  (75X)  of  the  arrivals.  (Obviously  75X  is  a  variable;  probably 
anywhere  between  60Z  and  90Z  is  alright.)  If  so,  then  all  correct 
arrivals  in  that  AA  event  may  be  added  to  the  total  correct  arrivals 
count;  and  all  incorrectly  associated  arrivals  should  be  subtracted.  T 
propose  that  a  correctly  associated  arrival  should  count  even  if  the 
phase  name  is  incorrect.  (Or  perhaps  it  should  get  1/2  or  0  weight  if 
incorrectly  identified.)  In  this  way  all  "perfect  bulletin"  events  are 
considered.  Any  remaining  arrivals  which  have  been  associated  into 
events  (which  are,  presumably  "false"  or  "bad"  events  and  which  will 
give  an  analyst  severe  trouble  to  clean  up)  will  be  subtracted  from  the 
total  correct  arrivals  count.  The  Association  Percentage  is  calculated 
as  (total  correct  arrivals  count)/(base  set  of  arrivals). 

This  measure: 

Doubly  penalizes  the  creation  of  false  events  and  incorrect 
association  to  good  events. 

Rewards  correct  identification  of  later  phases. 

Does  not  depend  on  location  error  which  is  highly  variable  depending 
on  event  size  and  travel-time  residuals. 

Can  be  easily  calculated  automatically. 

Perhaps  too  severely  penalizes  splitting  large  events,  especially  an 
equal  split. 
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INTRODUCTION 


The  problem  of  verifying  chat  magnitude  estimators  for  events  which 
have  been  created  by  an  automatic  association  (AA)  routine  are  based  on 
a  reasonable  subset  of  detecting  and  non-detecting  stations  has  been 
considered  by  Elvers  (1980).  She  developed  a  "plausibility  function" 
for  the  amplitude  distributions  which  is  essentially  proportional  to  the 
likelihood  of  observing  a  particular  configuration  of  amplitudes  when 
the  event  is  assumed  to  have  occurred. 

In  this  discussion,  we  attempt  to  improve  several  aspects  of  this 
procedure.  As  a  computational  improvement  we  also  have  developed  a 
Newton-Raphson  procedure  for  computing  the  maximum  likelihood  magnitude, 
assuming  the  signal  and  noise  variances  are  known.  This  procedure  runs 
much  faster  than  the  search  techniques  used  by  previous  authors. 

We  also  give  a  method  for  discarding  single  stations  with  outlying 
amplitudes  (I,  see  Figure  1)  based  on  the  "influence"  they  exert  on  the 
magnitude  estimator.  This  replaces  the  Elvers  procedure  of  comparing 
the  separate  components  of  the  likelihood  for  the  purpose  of  determining 
the  influential  stations.  Our  procedure  for  discarding  stations  is  also 
apparently  not  so  likely  to  "destroy"  the  event  in  the  case  that  a  large 
number  of  stations  are  "down"  and  fail  to  report,  but  do  not  inform  the 
bulletin  analysis  center  of  this  fact. 

Slunga  (personal  communication,  1982)  has  pointed  out  that  the 
"dynamic"  amplitude  criteria  of  Elvers  seem  to  erroneously  discard 
arrivals  from  an  event  on  the  basis  of  amplitude  anomalies.  On 
occasion,  although  the  amplitudes  are  indeed  anomalous,  the  probability 
that  the  arrival  time  could  agree  by  accident  seems  even  smaller,  so 
that  it  seems  unreasonable  to  totally  discard  the  arrival  from  the 
event . 

To  help  remedy  this  paradox  we  propose  to  measure  the  overall  event 
plausibility  with  a  "kinematic"  likelihood  ratio  (II)  in  parallel  with 
the  dynamic  one.  If  the  dynamic  criteria  suggest  that  an  observed 
amplitude  is  unreliable  but  the  kinematic  criteria  are  favorable,  then 
we  suggest  that  the  appropriate  procedure  in  most  cases  is  to  assume 


that  there  is  some  error  in  the  amplitude  measurement  and  to  recalculat 
the  magnitude,  after  making  an  appropriate  change  in  the  amplitude  data 
but  to  continue  to  use  the  reported  arrival  time  for  estimation  of 
location. 

Finally,  we  have  developed  a  procedure  for  tying  these  new 
computational  techniques  together.  In  general  terms  this  procedure  is 
not  derived  from  theoretical  considerations  but  embodies  common  sense 
ideas  about  the  causes  of  erroneous  amplitude  measurements. 

It  is  important  to  emphasize  that  these  techniques  cannot  simply  be 
inserted  into  an  automatic  association  program  as  mathematical  routines 
like,  for  example,  the  cosine  function.  Instead  the  routines  contain 
parameters  which  are  specific  to  the  network,  detectors,  and  methods  of 
analysis.  These  parameters  must  be  determined  by  careful  statistical 
analysis  of  network  bulletins  which  are  relatively  error-free,  and  by 
analysis  of  false  events  produced  by  the  automatic  association  program 
in  which  they  are  to  be  used. 
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FUNCTIONAL  DESCRIPTION 


In  Figure  1  we  see  a  flow  chart  of  the  proposed  procedure.  The 
Roman  numerals  I-Il  denote  the  computational  techniques  referred  to  in 
the  Introduction. 

At  the  top  of  Figure  1  we  see  that  the  event  is  submitted  to  the 
program  for  analysis.  Obviously  only  those  events  should  be  accepted 
for  which  the  program  is  prepared.  For  example  we  may  anticipate  that 
the  first  version  of  the  program  will  be  prepared  to  work  only  on 
initial  P  waves.  If  the  event  is  made  up  of  Pg ,  S,  Lg ,  etc.  then  the 
algorithm  has  nothing  to  contribute.  On  the  other  hand  an  event  with 
only  one  P  wave  and  without  an  amplitude  reading  may  be  submitted.  The 
single  detection  might  be  unlikely  because  many  stations  with  lower 
thresholds  nearby  the  detecting  station  did  not  detect.  When  this  P 
wave  is  discarded  then  the  event  might  be  unreliable  on  kinematic 
grounds . 

The  event  data  is  then  submitted  to  a  subroutine  (I)  which  computes 
the  Ringdal  magnitude  and  detects  outlying  magnitudes.  In  Appendix  I  is 
a  discussion  and  listing  of  a  subroutine  (amp)  which  rapidly  computes 
the  Ringdal  magnitude  given  station  amplitudes  and  thresholds  in  terms 
of  magnitudes,  and  assuming  that  the  signal  and  noise  variances  are 
known.  Before  this  subroutine  can  be  used,  of  course,  distance- 
amplitude  relations  must  be  used  to  transform  amplitudes  and  thresholds 
t o  magn i t ud e s . 

The  amplitude  patterns  can  be  subdivided  roughly  into  n^  observed 
amplitudes,  Oj  above  the  noise  threshold  and  with  arrival  times  but  not 
amplitudes  reported,  and  n^  below  the  noise  threshold  and  not  detected. 
Of  course  other  subcategories  such  as  clipping  might  be  added.  For  the 
general  unknown  magnitude  case  the  log- 1  ike  1 ihood  as  given  by  Elvers 
(1980)  would  be  of  the  form: 
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is  a  standard  value  depending  on  the  mean  noise  level  D.  with  variance 

2  J  2 
o-  ,  the  theoretical  magnitude  m  with  observed  variance  o  ,  and  a 
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signal-to-noise  constant  C.  Since  m  is  the  only  unknown,  (1)  can  be 

maximized  using  the  Newton-Raphson  or  scoring  algorithms,  which  leads  to 
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the  maximum  likelihood  estimator  m  and  its  estimated  variance  0-  (see 

m 

Appendix  I).  .The  above  approach  is  essentially  that  followed  by  Elvers 
(1980). 

As  a  modification  to  the  above  Elvers  (1980)  suggests  comparing  the 
value  of  the  maximized  log  likelihood,  say  In  L(m)  with  some  arbitrary 
threshold,  with  the  first  component  of  (1)  modified  to  behave  like  an 
estimated  interval  probability.  If  the  total  plausibility,  say  In  L(m) 
is  too  small,  then  the  event  can  be  rejected.  The  other  decision  that 
can  be  made  at  this  point  in  the  Elvers  procedure  is  to  discard  a  single 
station  value  based  on  a  single  station  component  of  In  L(m).  For 
example  the  decision  to  reject  a  station  j  which  does  not  observe  is 
based  on  the  value  of  ln($(-Zj)),  which  is  the  probability  of  observing 
a  value  below  the  noise  threshold  at  the  jth  station. 

The  problem,  mentioned  by  Elvers  (1980),  with  this  approach  is  that 
the  total  plausibility  depends  in  part  upon  point  values  of  the  density 
and  in  part  upon  interval  probabilities  which  are  integrated  densities. 
Thus  it  is  likely  that  eliminating  an  event  purely  on  this  criterion 
will  unnecessarily  eliminate  events.  The  problem  is  magnified  when  one 
uses  the  Elvers  (1980)  procedure  to  eliminate  single  stations;  purely  on 
the  basis  of  their  plausibility  values. 

As  a  better  method  for  detecting  station  outliers,  we  suggest 
recalculating  the  magnitude  estimator  with  each  station  value  missing  to 
obtain,  say  for  i“l,...,n,  where  m^  ^  denotes  the  magnitude 

estimated  when  station  i  is  left  out  of  the  calculation.  These  can  be 
compared  with  the  original  estimator,  say  m,  to  determine  the  "influence 


of  station  i. 


For  example,  a  common  measure  of  influence  is  "Cook's 
Hist  mice",  defined  in  this  case  by: 


t-O 


A  A 


(2) 


An  heuristic  procedure  suggested  by  Cook  is  to  reject  station  i  when 
z(_i)  mt>ve8  further  away  from  zero  than  z(a/2),  where  <*  is  some 
arbitrary  probability  value. 


As  an  example  of  how  the  two  procedures  for  amplitude  verification 

might  work  in  practice,  consider  the  11  station  "event",  shown 

schematically  in  Table  1.  We  note  first  that  taking  the  simple  mean  of 

the  observed  magnitudes  yields  "m  -  4.04,  whereas  the  maximum  likelihood 

estimator  ism"  3.78  with  estimated  standard  deviation  o*  *  .12  so 

m 

that  the  censored  observations  definitely  pull  the  estimator  down. 


In  order  to  determine  whether  any  of  the  station  values  m.  are 
outliers,  we  may  use  the  analogue  of  Cook's  distance  given  in  equation 
(2).  The  estimated  deleted  station  magnitudes  m^_.j  for  i"l,2,...,ll 
are  shown  in  Table  1  and  we  note  that  the  greatest  change  is  produced  by 
omitting  station  6(m_g  ■  3.89)  which  is  the  value  known  to  be  below  its 
noise  threshold  »  3.  The  elimination  of  this  station  then  moves  the 
maximum  likelihood  estimator  back  towards  the  mean. 
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TREATMENT  OF  STATIONS  WITH  ANOMALOUS  AMPLITUDES 

Returning  to  Figure  1  we  must  ask  what  course  of  action  to  take  when 
an  anomalous  station  amplitude  is  discovered  by  procedure  I.  First  we 
determine  if  the  most  anomalous  station  has  reported  a  detection  or  if, 
as  in  Table  1,  it  has  simply  failed  to  report  at  all,  implying  that  the 
signal  level  is  below  the  noise  level.  In  the  latter  case  we  follow  the 
approach  suggested  by  Elvers  (1980)  and  simply  discard  the  station. 
However,  we  may  look  ahead  a  bit  and  note  that,  unlike  F.lvers  we 


TABLE  1 


MAXIMUM  LIKELIHOOD  ESTIMATORS  FOR  MAGNITUDE 


Case  I 

i 

Data 

A 

“-i 

A 

0» 

m 

InL 

*(-i> 

Station 

1 

-  B 

3.74 

.13 

1.03 

.32 

D 

3.80 

.12 

.96 

-.25 

B 

3.67 

.13 

2.27 

.89 

B 

3.74 

.13 

1.03 

.32 

5 

3.70 

.13 

1.51 

.60 

* 

6 

B> 

3.89 

.17 

6.14 

-1.00 

m 

3.80 

.12 

3.06 

-.23 

4.2 

3.79 

.12 

2.87 

-.14 

B 

3.9 

3.81 

.12 

3.20 

-.29 

4.5 

3.78 

.12 

2.73 

-.05 

9 

5.0 

3.78 

.12 

2.68 

-.00 

immediately  test  for  the  kinematic  probability  of  this  event  (II)  and 
if  the  probability  of  the  event  being  real  is  less  than  say  0.95  then  we 
exit  the  routine  with  a  bad  event  flag.  That  is,  non-detection  by  this 
station  suggests  that  there  is  something  wrong  with  the  event,  and 
unless  the  kinematic  probability  is  high  the  event  is  discarded.  On  the 
other  hand,  if  the  kinematic  probability  of  the  event  is  high,  then 
probably  there  was  something  wrong  at  the  station,  and  we  return  to 
further  analyze  the  event  without  considering  this  station.  Since  in 
practice  many  stations  fail  to  report  due  to  operational  consideration 
we  should  probably  set  a  less  restictive  threshold  for  the  non-detecting 
stations,  that  is,  we  should  tend  to  reject  them  first.  For  example,  if 
for  a  non-detecting  station  z^  ^  -0.7  reject  it  (see  Table  I)  but 

otherwise  require  z^_^j  1.0  for  rejection.  Probably  we  should  have 

no  anomalous  non-detecting  station  remaining  before  considering  other 
stations.  Further  work,  perhaps  empirical  is  needed  to  best  determine 
those  thresholds. 

If  the  anomalous  station  did  detect  we  ask  if  an  amplitude  was 
measured  or  not.  If  no  amplitude  was  measured  then  we  must  have  the 
situation  where  this  station  could  not  have  been  expected  to  detect; 
typically  a  station  in  a  shadow  zone  to  a  small  event.  Again  we  discard 
the  station  and  check  the  kinematic  probability.  Since  we  have  by  this 
time  eliminated  many  of  the  incorrectly  non-reporting  stations,  the 
event  magnitude  should  not  be  biased  too  low  so  that  stations  should  not 
be  incorrectly  thrown  out  at  this  point. 

Finally  we  have  the  case  where  there  is  an  amplitude  measurement  at 
the  anomalous  station:  either  it  is  too  large  or  too  small.  In  either 
case  we  assume  that  some  blunder  has  been  made  and  that  if  a  true  signal 
has  been  detected  then  the  amplitude  has  been  recorded  or  transmitted 
incorrectly.  Thus  we  simply  change  the  amplitude  measurement  (not  in 
the  original  files  of  course)  to  state  that  the  signal  was  greater  than 
the  noise.  Again  we  proceed  through  the  kinematic  criterion.  Note  that 
if  this  event  comes  around  again  then  this  detection  will  not  have  a 
measured  amplitude  so  that  if  it  is  unlikely  that  this  station  could 
have  detected  the  event  at  all  the  arrival  will  be  discarded. 


KINEMATIC  REJECTION  CRITERIA  (II) 


The  kinematic  criteria  for  the  event  existence  are  to  be  established 
as  follows.  First,  it  is  necessary  to  establish  event  types  which  will 
consist  of  the  number  of  array  stations  and  the  number  of  non-array 
stations  detecting.  These  stations  may  also  have  to  be  broken  down  by 
the  criteria  of  analyst  and  automatic  detection.  For  each  type  of  event 
we  may  determine  the  ratio  of  good  events  to  false  events  of  each  type 
in  each  day.  This  ratio  then  gives  directly  the  probability  that  an 
event  of  each  type  is  real. 


PROCESSING  AFTER  DYNAMIC  CHECKING 

If  any  detections  have  been  discarded  then  the  reduced  set  of 
arrivals  should  be  used  to  generate  a  new  trial  epicenter.  If  the  same 
final  set  of  arrivals  should  result  then  the  iteration  should  be 
suspended.  Care  needs  to  be  taken  to  avoid  an  infinite  loop. 
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APPENDIX  I 


We  have  developed  a  Newton-Raphaon  algorithm  for  maximizing  In  L(») 
by  noting  that: 
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where  z.  given  by  the  equation  in  the  main  text  with: 
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memorandum 


J.  Goncz 


FROM: 


R.  Bland  ford  v!  & 


SUBJECT:  A  Better  Approach  to  Handling  Amplitude  Data  in  AA 

DATE:  19  April  1983 

1.  Compute  maximum  likelihood  m.  using  the  observing  stations  plus  a 
fixed  set  of  about  20  reliable  stations.  An  m^  is  calculated  with 
each  station  omitted  in  turn.  Assume  an  a-priori  o  ■  0.35  m^ 
throughout. 

2.  Analyze  each  of  the  above  stations  in  screen  as  we  do  now,  get  the 
total  number  of  points,  then: 

i.  If  the  station  did  not  detect  but 

p(det)  >  .95,  -1  point 
p(det)  >  .99,  -2  points 

Comments:  If  p(det)  >  .95  at  many  stations  most  stations  will 

have  detected  if  the  event  is  big  and  loss  of  points  for  the 
occasional  random  non-detection  will  not  matter.  If  the  event 
is  weak,  only  a  few  close  stations  will  have  p(det)  >  0.95  so 
the  occasional  random  non-detection  is  again  unlikely  to  lose 
points . 

ii.  If  the  station  did  detect  and  if,  based  on  noise  and  m^  alone 

p(det)  <  .01, -(all  points  due  to  station) 

Comments:  If  the  event  is  big  there  will  be  so  many  points  that 

occasionally  losing  one  will  not  be  serious.  If  it  is  a  small 
event  and  just  trapped  (correctly)  the  one  distant  station,  then 
the  probability  of  that  happening  for  any  particular  one  out  of 
the  20  distant  stations  may  be  smaller  than  .05  so  for  that 
reason  we  set  this  threshold  at  .01  so  that  only  about  1/5  of 
these  events»will  be  incorrectly  rejected.  It  is  important  to 
reject  all  the  points  due  to  the  station,  i.e.  those  due  to 
azimuth  and  slowness  and  those  due  to  local  flags  and  associated 
S  phases. 

iii.  If  the  station  did  detect,  and  reports  an  amplitude  Aq,  and  if, 
based  on  noise  and  m^  alone, 

p(det)  >  .01  and 
p(amp  >  Ao)  <  .01  or 
p(amp  <  Aq)  <  .01  then 


M 
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discard  amplitude,  A  (are  left  with  estimated  noise)  and  return 
to  (1)  0 

Comments:  That  is,  compute  all  the  negative  points  and  m^ 

values  over  because  the  bad  amplitude  probably  distorted  all  the 
m^  values.  Note  that  this  procedure  does  not  "fiddle"  with  an 
event.  It  either  results  in  complete  acceptance  or  discards  the 
event  so  that  the  kinematic  processing  has  a  new  chance  to  do 
better. 
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